The Weyl Series and the Trace formula: Can we add them? 
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Periodic orbit expressions for the density of states lead to spurious results when directly used to 
calculate quantities of thermodynamic interest. This is because the trace formula is usually valid 
only for large energies while the calculations make use of it at all energies. We present a prescription 
for circumventing this problem by isolating contributions which arise from the inaccurate low-energy 
behavior, and are spurious, from other subdominant contributions which are physical and are not 
contained in the Weyl series. The method is tested by analyzing the fermionic grand partition 
function for a disk billiard and may readily be extended to other dynamical systems. 
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In many-body systems such as atomic nuclei, metal 
clusters, quantum dots and traps with a dilute gas of 
bosons, much of the physics may be unraveled from the 
quantum description of the single particle motion in an 
appropriate mean field. Much effort has gone into de- 
veloping intuitive, semiclassical descriptions of single- 
particle quantum dynamics One ubiquitous con- 
cept is the division of the density of states into two com- 
ponents, each with a distinct semiclassical interpretation. 
The first varies smoothly with energy and other param- 
eters and is related to the geometrical properties of the 
classical phase space; for billiard systems this is the well- 
known Weyl series. The second component is an oscilla- 
tory function and is related to the dynamics of the clas- 
sical system through various trace formulas. Balian and 
Bloch H showed that the decomposition can be made 
exact at the price of expressing the result in terms of 
certain integral representations. In practice, this is not 
very helpful and one is led to the expedient of expanding 
these out as asymptotic series. This leads to the Weyl 
term plus an infinite sum over periodic orbits, each one 
of which has its own asymptotic expansion in powers of 
H. (Explicit higher order terms in the trace formulas for 
more than one dimension have been worked out only in 
a few trial systems [||6|.) We will refer to the sum over 
the periodic orbits to leading order in % as the trace for- 
mula. In the resulting approximation, however, it is not 
immediately clear to what extent this leads to spurious 
results in integrated quantities derived from the trace for- 
mula. The focus of this paper is to investigate this ques- 
tion, given a representation of the density of states as an 
infinite asymptotic Weyl series plus the trace formula. 
Certainly, this is the form in which the semiclassical re- 
sults are most commonly expressed and used We 
would like to know to what extent we can meaningfully 
use this as it stands: for example in using it to calcu- 
late quantities of thermodynamic interest by integrating 



over it with various weight functions. The short answer 
is that we can not. The Weyl formula will typically have 
terms which cannot be integrated directly from zero en- 
ergy. This problem is, however, easily circumvented by 
judicious integrations by parts [Q. A more serious prob- 
lem, investigated here, is related to the fact that the trace 
formula is usually inaccurate at small energies and can 
lead to grossly incorrect results for the integrated quan- 
tities. We will show that this problem can be overcome 
by identifying a spurious component, which is due to the 
incorrect low-energy behavior. There is another compo- 
nent that is subdominant but contains the real, physical 
information about the periodic orbits. 

In this paper we focus on the disk billiard, in part be- 
cause its Weyl expansion is known to very high order || 
and also because of its physical importance, for exam- 
ple in quantum dots §. However, at the end we mention 
analogous results for a variety of systems in both two and 
three dimensions. For the disk, it is shown in || that the 
divergence of the higher-order terms in the Weyl series is 
controlled asymptotically by the shortest accessible peri- 
odic orbit through the phenomenon of resurgence JlO|,[llJ • 
The intricate link between the trace formula and the 
Weyl series has also been studied by Cartier and Voros 
fll2| , [l3| |. Another interesting system to study in this re- 
gard is the harmonic oscillator JjJ where there is an exact 
trace formula. Since it can be found by resummation of 
the Weyl series, we can work with either form. There are 
no spurious terms although there are still contributions 
from the periodic orbits which are exponentially small in 
temperature, as we will find for the disk. 

We begin with the following decomposition of the 
quantum density of states: p(E) = J2n$(E ~ E n ) w 
~p(E) + Sp(E). The Weyl part ~p{E) is given in terms 
of reciprocal powers of E as well as delta functions and 
their derivatives. (The corresponding expression for the 
partition function is simpler and is presented below.) 



The trace formula for the oscillating part 5p(E) (with 
h 2 /2m = 1 and disk radius R — 1) is given by @jH 

oo [v/2] ^ 

5p{E) « E E sin(VSi vw + 9 V ), (1) 

v=2 w=l 

valid for E 3> 1. Each pair of integers (v, w) represents 
a one-parameter family of orbits, related by continuous 
rotations, where v is the number of vertices (i.e., reflec- 
tions from the boundary of the disk) and w is the wind- 
ing number. L vw — 2vsin(-7rw/v) is the length of an 
orbit, A vw = /vw^vw /\/87rv 2 its amplitude and 9 V = 
(3?r/4 — 3v7r/2) its phase. The degeneracy factor / vw is 
unity for v = 2w, otherwise it equals 2. We now consider 
using this expression to calculate physical observables of 
interest. As an example we consider the fermionic grand 
partition function [ p"5| for non-interacting particles, from 
which all thermodynamic quantities can be derived: 

/>oo 

lnZ G (a,0) = p{E) In [1 + exp(a - 0E)] dE . (2) 
Jo 

Here is the reciprocal temperature and a/0 = p, is the 
chemical potential. However the considerations of this 
paper will apply to any quantity expressible as an integral 
over the density of states, such as the total energy of a 
many-fermion system or the integrated density of states. 

To proceed it is convenient to introduce the one-body 
partition function 

/>oo 

Z(s)= p{E)exp{-sE) dE . (3) 
Jo 

The parameter s can be thought of as an inverse tempera- 
ture for the one-body problem. In what follows we will al- 
low it to be complex, and for the present purpose it is just 
used as a mathematical variable. Noting that the two- 
sided Laplace transform of In [1 + expja)] is 7r/ssin(7rs) 
and using the convolution theorem |p_ T|| we can rewrite 
Eq. (^) in the form 

hiZ G (a,0) =±-( dse as { % , Z(0s)\ . (4) 
^ l Jc I s sm(7rs) J 

The contour C runs from — ioo to ioo with the real part 
of s between and 1, although obviously this can be de- 
formed within the rules of complex calculus. The study 
of the grand partition function has the nice feature of 
being directly related to physical observables. It has also 
the following mathematical advantage. It is well known 
that the trace formula for the single-particle density of 
states is generally divergent. However the finite temper- 
ature mitigates this. The number of orbits in the disk 
increases as a power law with orbit length while the tem- 
perature causes an exponential suppression with temper- 
ature. The net effect is a convergent sum. 
Substituting ~p(E) into (||), one obtains 



s ^ T[r 2) 

with Co = 1/4, ci = — 7r/4, ci = 1/6, C3 = 7r/256, etc. 
Berry and Howls || tabulate the coefficients up to C31. 
This expression for Z may be substituted in Eq.(^) to 
obtain the Weyl series for the grand potential. We find 

h[Z^(a, (3) = -| / dE\n(l + e a - E ) (6) 
P Jo 

N J"-— 1 

y c 2n+l on-1/2 d / d£ :_Lwi + e «-S). 
^T(n+l/2f da" J ^E 1 ' 

This series is asymptotic and diverges if taken to infinite 
order; it may therefore only be summed up to a maximum 
value N of the index n. This is shown in Fig. 1, where 
the difference 

A(a,/3) = In Z G (a,f3)-\^ZG~(a,0) (7) 

between the quantum and the Weyl value is plotted ver- 
sus N (circles) for selected fixed values of a and 0. 
For /i = a/0 ^> 1, A settles down to a plateau after 
a few terms of the series, giving the optimum result. 
The asymptotic nature of the series (^) is apparent for 
a/0 = 1 in the bottom of Fig. 1, where |A| is seen to 
increase rapidly with the inclusion of higher order terms. 
(The same happens for the other cases, but at higher val- 
ues of AT.) We see from this figure that the series yields 
a plateau at a nonzero value, indicating that (^|) devi- 
ates significantly from the true quantum result. This is 
reminiscent of the phenomenon discussed by Balian et 
al. Jl6| where it is shown how the error in using a given 
asymptotic expansion can be dominated not by the least 
term in the expansion used but rather in neglecting some 
weaker subdominant saddle contribution. Clearly what is 
missing here is the contribution from the periodic orbits. 

We proceed by substituting into (||) the oscillating den- 
sity of states as given by (|]). This gives a spurious result. 
For example, in the limit — we are left with simple 
Fresnel integrals. Upon doing the sum over orbits we 
get the finite contribution (In 2 - ?r 2 /24) ln(l + e Q )/4\/2. 
However, we know that in this limit the Weyl contri- 
bution is exact and clearly adding this amount would 
destroy the agreement. Therefore, we should not take 
the prescription of substituting the oscillating density of 
states into (|^) too literally. 

Rather, we proceed as before by first calculating the 
contribution to the one-body partition function. Chang- 
ing variables to k = \f~E~ we have the following integral 
to do for each orbit 

/•OO 

/ dkVksm(kL + 9)exp(-sk 2 ), (8) 
^0 



where we have temporarily suppressed the v and w in- 
dices. We decompose the sinusoid into two exponentials 
and then separately analyze the two integrals. There is 
an end-point contribution from the lower limit which can 
be determined from Watson's lemma There is also 
a saddle at either of k = ±iL/2s which may or may 
not contribute depending on the phase of s; there is a 
Stoke's phenomenon when 3s = 0. The final result valid 
for Isl < 1 is fH 



V8v 2 ^ 22»n! \Ll w 



S7< \~ ST f H V/2 ^(4n + l)!! ( s 



±i(6» v — tt/4) 



(9) 



(10) 



where the upper and lower signs in the phase of ( |Io|) hold 

for 3s < 0. When 3s = 0, one should take the mean of 
the two expressions [ p"(i| . The sum over n in (^J) is asymp- 
totic and must be truncated at some maximum value N s 
that depends on the orbit. The series ([)]) comes from the 
end-point analysis. When substituted into (Q), its first 
term is the unwanted contribution at (3 = mentioned 
earlier. In fact, the entire series is unwanted. It results in 
power series contributions to \tiZq, but we have already 
shown in Eq. (||) that its power series is correctly gener- 
ated by the Weyl series alone. Therefore this entire series 
is spurious — hence the subscript s — and is neglected 
in what follows. This will be justified more rigorously 
below. The second contribution (10), by contrast, comes 
from the saddle-point analysis. We will see that it does 
contain the physical information about the contributions 
of the periodic orbits — hence the subscript p. It is 
interesting to note that ( [To| ) is exponentially subdomi- 
nant compared to (^|) and yet carries all the relevant in- 
formation. Higher-order saddle-point corrections to JlO| ) 
have been worked out Jl9| , but there is no reason to in- 
clude them here because they correspond to the missing 
higher-order corrections to the trace formula (Jj]) which 
are unknown. 

The origin of the spurious scries (^|) is the low-energy 
region, where the trace formula is inaccurate. Removing 
this spurious series could be thought of as subtracting 
from (0) a series of the form Sp s (E) — J2n=o a n ^ n \E). 
Requiring that all moments of the corrected trace formula 
be zero fixes the coefficients a n and leads to subtraction 
of (^|) from SZ(s). Another argument for ignoring SZ s (s) 
is as follows. The derivation of ([!]) assumes energy is 
large. We are therefore free to replace this expression by 
any other which is asymptotically the same. However, 
the expansion (^) is not invariant under such a change 
(unlike (|l(J)) and we therefore conclude that it cannot 
contain any meaningful information. It is even possible 
to find a version of Sp(E) whose spurious contribution 
is identically zero by inverse Laplace transforming the 



physical expression (|l0|). This gives Bessel functions |7j 
whose asymptotic expansions are precisely Sp(E). How- 
ever for small energy they behave more smoothly and, 
most importantly, do not lead to spurious structure. We 
interpret the result as a regularised trace formula. 

We now calculate the contribution of the periodic or- 
bits to the fermionic grand partition function by sub- 
stituting SZp(s) into (^) (while omitting SZ S ) and do- 
ing the contour integral numerically. We call the result 
(S\nZc) p and show it in Table 1 for a few cases together 
with A(a,/3). Clearly this analysis captures the part 
missing from the Weyl series. We stress that this would 
certainly not be the case had we included the spurious 
series. We also show this data in Figure 1 as the dashed 
horizontal line. For comparison we also show the result 
of substituting (|l|) directly into (||), doing the integrals 
numerically and subtracting off the spurious series found 
from substituting (||) into the relation (ji|). The results 
are shown as crosses where the horizontal axis represents 
N s in ([)]). The asymptotic result is well captured by 
(SlnZ G ) p . 

Thus, the contribution of periodic orbits to the grand 
potential accounts for the difference between the quan- 
tum and the Weyl calculations in a consistent manner, 
but only after discarding the spurious contributions. It 
is interesting to note that the periodic orbits represent 
a much more important contribution to In Zq than to Z 
itself, for reasons we do not have space to discuss here 
0. The integral (||) has saddles at s = ±iL/2y / a//3 and 
the resulting stationary phase analysis leads to a known 
form |^0| valid in the limit of large a. However, by do- 
ing the integral exactly we can relax this constraint and 
therefore have a more general result. (Both approaches 
fail, however, when p becomes much smaller than one.) 

This sort of analysis can immediately be generalized 
to other classes of periodic orbits. For example, a 
typical billiard orbit contributes to trace formulas as 
AE n l^ sm(v r EL+air/2+mr/4). Examples are diffractive 
orbits in two dimensions, isolated orbits, the disk orbits 
considered here, diameter orbits in the sphere and polyg- 
onal orbits in the sphere for which n = (—3,— 2,— 1,0,1) 
respectively. In its contribution to Z(s), we can identify 
a component which comes from the end point and is spu- 
rious and another which comes from a saddle point and 
is physical. The latter is approximately 



»/2+l p -L 2 /is 



„(n+3)/2 



D ±i(a+n)ir/2 



(11) 



which agrees with (|l(]) for the case n = — 1 with a = 
2 — 3v and substituting in the amplitudes from ([!]) . This 
can also be extended to potential systems, scaling or oth- 
erwise, although the former is somewhat simpler. 

In summary, a careful asymptotic analysis of integrals 
over approximate trace formulas reveals that the physi- 
cally relevant contributions from the periodic orbits are 



contained in exponentially subdominant terms, whereas 
the dominant terms that come from the inaccurate low- 
energy behavior of Sp(E) are spurious and must be dis- 
carded. Adding only the subdominant terms to the Weyl 
series gives good numerical agreement with quantum- 
mechanical results. 
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a 


p 




A 




10 


1.0 


4.2499 


0.2880 


0.2549 


10 


0.5 


12.6489 


0.2144 


0.2100 


10 


0.1 


97.0390 


-0.0075 


-0.0071 


1 


1.0 


0.0083 


0.0148 


0.0152 



TABLE I. The logarithm of the grand partition function, 
the difference from the Weyl approximation and the contribu- 
tion from the periodic orbits (with the spurious part removed) 
for a selection of values of a and (3. 
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FIG. 1. The logarithm (base 10) of A versus the order in 
the asymptotic expansion (circles). The dashed line is the 
periodic orbit prediction (Sin Zc) p - The crosses result from 
directly using the expression (0) in the energy integral (^) and 
then explicitly subtracting off the spurious component versus 
N s in 



